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Abstract 

We study the dynamics of charged pions in the nuclear medium via the 
ratio of differential 7r“- and Tr^-spectra in a coupled-channel BUU (CBUU) 
approach. The relative energy shift of the charged pions is found to correlate 
with the pion freeze-out time in nucleus-nuclens collisions as well as with the 
impact parameter of the heavy-ion reaction. Fnrthermore, the long-range 
Coulomb force provides a ’clock’ for the expansion of the hot nuclear system. 
Detailed comparisons with experimental data for Au + Au at 1 GeV/A and 
Ni + Ni at 2.0 GeV/A are presented. 
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1 Introduction 


The meson-baryon dynamics in relativistic nnclens-nnclens collisions is of particnlar 
interest since one expects to learn abont the hadron properties and abnndancies in 
hot and dense nnclear matter. Whereas dileptons from decays or 7r'''7r“ annihi¬ 
lation provide information on the timescale of 1.3 fm/c in the high density phase 
B B B B B0’ pion correlations can be used to learn about the freeze-out volume 
of the pions during the expansion phase. Whereas the latter method has been ex¬ 
tensively explored in nucleus-nucleus collisions at AGS and SPS energies [Q, pion 
correlation measurements are less promising at SIS energies due to the much lower 
pion multiplicities achieved at 1 - 2 GeV/A. Alternatively, one might here use the 
kinematic shift of 7t~ to tt’*' spectra due to the long range Goulomb force p| in order 
to learn about the baryon expansion dynamics at the production time of those pions 
that no longer interact with the baryons and are identified experimentally [y, p!0[| . 

In this paper we investigate more systematically the pion dynamics in the 
time dependent Goulomb field generated by all charged hadrons in nucleus-nucleus 
collisions from 1-2 GeV/A for systems at SIS. Our analysis is carried out within the 
coupled-channel BUU (GBBU) transport approach which is briefly described in 
Section 2. In Section 3 we analyse the correlation between Goulomb effects and the 
pion production time as a function of the pion momentum in the cms and study the 
sensitivity of the ti~ spectra to the reaction geometry of the heavy-ion collision. 
Detailed comparisons with experimental data for the ti~ ratio are, furthermore, 
presented for Au + Au at 1 GeV/A and Ni -|- Ni at 2.0 GeV/A. Section 4, hnally, 
concludes our study with a summary and a discussion of open problems. 


2 The CBUU-Model 


2.1 Basic equations 


For our present study we use the GBUU-transport-model [|^ to describe the time 
evolution of relativistic heavy-ion collisions. In this approach, apart from the nucleon 
and the A(1232), all nucleon resonances up to masses of 1.95 GeV/c^ are taken 
into account as well as the mesons tt, rj, p and a. While the vr-, rj- and p-mesons 
correspond to physical particles, the cr-meson is introduced to describe correlated 
pion-pairs with total spin J = 0. For the baryons as well as for the mesons all isospin 
degrees of freedom are treated explicitly. The hadrons included in our model obey 
a set of coupled transport equations for their one-body phase-space distributions 

mpA) [|T10,[1,[1]: 
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The l.h.s. of eq. (|^ represents the relativistic Vlasov-equation for hadrons with 
charge qi moving in a scalar momentnm-dependent held Ui{f,pi) H as well as 
in the electromagnetic potential Vc'(T), where r and pi stand for the spatial and 
momentnm coordinates of the hadrons, respectively. The effective mass m*(r, pi) in 
eq. (0) inclndes the restmass of hadron i as well as a scalar mean-held potential 
Uiif, pi), i.e. 

m*ir, pi) = +Uiir, pi). (2) 


Since in our model the mesons will be propagated as free particles with respect 
to the nuclear interaction, their ehective masses are equal to their restmasses, i.e. 
Ui{f,pi) = 0 for mesons. The r.h.s. (i.e. the collision integral) of eq. (D describes 
the changes of /j(r,due to two-body collisions among the hadrons (hi + h 2 <->■ 
hs -|- hf) and two-body decays of baryonic and mesonic resonances. Furthermore, 
is the in-medium collision rate, /i = 1 — /i (* = are the Pauli- 
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blocking factors for fermions and ni 2 is the relative velocity between the colliding 
hadrons hi and ^2 in their center-of-mass system. In the collision integrals describing 
two-body decays of resonances one has to replace the product (relative velocity x 
cross-section x / 2 ) by the corresponding decay rate and to introduce the proper 
fermion blocking factors in the hnal channel. The factor g in (^) stands for the 
spin degeneracy of the particles participating in the collision whereas ^ 2 , 3,4 stands 
for the sum over the isospin degrees of freedom of particles 2, 3 and 4. We include 
the following elastic and inelastic baryon-baryon collisions, inelastic meson-baryon 
collisions as well as meson-meson collisions (for a detailed discussion of the cross 
sections employed cf. 
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where R and R' represent all baryonic resonances up to a mass of 1.95 GeV/c^ |n 
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2.2 The Coulomb Potential 


Charged baryons and mesons are propagated in the electromagnetic field generated 
by the 4-current of all charged particles. Since our numerical studies did not indicate 
any sensitivity to the hadron trajectories from the time dependent magnetic fields, 
we only discuss the effects of the time-like component (i.e. Vc) in the following. 
The Coulomb potential in our transport model is evaluated by solving the static 
Poisson-equation in three dimensions 

Vc{x, y, z) = Airahcpcix, y, z) (4) 

{a = 1/137) in every time-step of the calculation where all charged particles present 
contribute to the charge density pc{x,y,z). The Poisson-equation (||) is integrated 
by means of the Alternating-Direction Implicit Iterative (ADI-)algorithm OB. 
where the boundary conditions are obtained by a multipole expansion up to the 
quadrupole moment of the Coulomb potential. Since the mesons are propagated as 
free particles with respect to the nuclear mean-field, the Coulomb force 

Fc{r) = -qVrV,{r) (5) 

(besides the Lorentz force from the magnetic field) is the only force acting on a 
meson with charge q. As noted before, the curl of the spatial components of the 
vector potential is found to be small in the overlap regime of the colliding ions and 
can be neglected within the numerical accuracy achieved. For charged baryons F/ 
represents a force in addition to the nuclear force from the scalar mean-field potential 
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however, is of minor importance here. 


( 6 ) 


3 Coulomb analysis of pion spectra 

In the following we investigate the Coulomb-effects on charged pion spectra in heavy- 
ion collisions in the SIS energy range. We especially look at the correlation between 
the Coulomb-forces acting on charged pions and their freeze-out times and freeze- 
out densities, respectively. Finally, we present pion spectra and 'n~ j'n^ ratios from 
our transport calculation for Au-\-Au at 1 GeV/A and Ni-\-Ni at 2.0 GeV/A which 
are studied experimentally at SIS. 
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Figure 1: Number of particles in the transport calculation as a function of time for 
the reaction Au + Au at 1.0 GeV/A and impact parameter 6 = 0 fm: solid line: 
A(1232); dashed line: sum of all higher resonances; dotted line: pions; dash-dotted 
line: sum of baryonic resonances and pions. The short-dashed line shows the baryon 
density in a volume of 1/m^ in the center of the CMS of the colliding heavy-ions. 

3.1 Correlation between Coulomb-effects and pion produc¬ 
tion times 

In the high density phase of a heavy-ion collision in the SIS-energy range baryonic 
resonances like the A(1232), A(1440), A(1535) etc. are excited in nucleon-nucleon 
collisions or via multi-step processes in pion-nucleon or resonance-nucleon collisions. 
These excited baryonic resonances then can be deexcited in collisions with nucleons 
(e.g. A(1232)A ^ NN) or decay into nucleons and pions. If these decays take place 
in the high density phase of the heavy-ion reaction the emerging pions are likely to 
be absorbed again, but if the decays occur in the expansion phase of the reaction 
or in low density regimes, i.e. at the surface of the reaction volume, the pions may 
escape from the nuclear medium (freeze-out). 

This line of argument is supported by £g. where we show the number of baryonic 
resonances and pions present in our transport calculation as a function of time for a 
central (6 = 0 fm) Au -|- Au collision at 1.0 GeV/A. The short-dashed line indicates 
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the baryon density in a volnme of 1 fm^ in the center of the CMS of the two 
colliding nnclei. The population of the baryonic resonances sets in as the density 
in the overlap zone of the colliding heavy-ions increases and reaches its maximum 
at the end of the high density phase at about 15 - 20 fm/c. This holds for the 
A(1232) (solid line) as well as for the higher resonances (dashed and dotted lines). 
For the pious, however, we obtain a different picture. The number of pions increases 
monotonically with time and reaches its maximum (asymptotic) value at the end of 
the reaction, which is about a factor 1.2 higher than the maximum number of A’s 
in the high density phase. This implies that most of the pions observed at the end 
of the heavy-ion reaction stem from resonances decaying in the expansion phase of 
the reaction. 

In order to discuss this effect more quantitatively we show in £g. the number 
of asymptotically observed pions as a function of their freeze-out times for CMS 
momenta of 0.2 (solid line), 0.6 (dashed line) and 1.0 GeV/c (dotted line) for a Au + 
Au collision at 1.0 GeV/A in central collisions (6 = 0 fm). The production of pions 
in all three kinematical regions yields about 10% of its maximum at f ~ 15 fm/c, 
reaches its maximum between 20 and 30 fm/c and then decreases with time. Since 
the production of pions with low momenta requires less energy than the creation of 
high momentum pions, there are quantitative differences in the production times of 
pions in the different kinematical regimes^. Pions with low momenta are already 
produced in the beginning of the overlap of both nuclei while the production of pions 
with high momenta starts only if there is enough energy density piled-up in order to 
produce them. For a similar reason the relative production rate of high momentum 
pions falls off more quickly with time than that of pions with low momenta. In 
general, the hard pions stem dominantly from the heavier baryonic resonances in 


Light resonances live on average longer than heavy resonances since their decay 
width is smaller than that of heavy resonances. This is particularly true for the 
A whose momentum-dependent width leads to a rather long lifetime for low mass 
A-resonances. A second reason is that light resonances can be populated not only 
in the compressed phase of a heavy-ion reaction but also in less energetic nucleon- 
nucleon collisions in the expansion phase of the reaction while the creation of heavy 
resonances - equivalent to the production of hard pions - needs more energetic two- 
body collisions which are only present in the high density phase of the heavy-ion 
collision. Due to this fact we observe that the maxima of the pion-distributions in 
£g. ^ are shifted to smaller times t as the pions become harder. 

This can be seen more quantitatively in £g. ^ where the average production time of 
pions produced in central events is shown as a function of their momentum in the 
GMS of the two colliding heavy-ions for the reaction Au -|- Au at 1.0 GeV/A. The 
squares indicate the result of the GBUU-calculation and the solid line corresponds 
to a linear £t to the calculated results. From this hgure we can infer that soft pions 
are emitted on average approximately 10 fm/c later than hard pions with momenta 
of about 1.0 GeV/c. During this time the baryonic matter expands before the 
emission of the soft pions. Thus the production density for hard pions is on average 


^The ’production’ time here is defined as the time when the A decays into a nucleon and a free 
pion that propagates out to the detector without any further strong interactions. 
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Figure 2: Number of emerging pious as a function of their final production time 
for Au + Au at 1.0 GeV/A taking into account only central events (6 = 0 fm) for 
different CMS pion momenta; solid line: Pcms = 0.2 GeV/c, dashed line: pcms = 0.6 
GeV/c, dotted line: pcms = 1.0 GeV/c. 
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Figure 3: Average production time of pious as a function of their CMS momentum 
for Au + Au at 1.0 GeV/A taking into account only central events (& = 0 fm); 
squares: results of the CBUU-calculation, solid line: linear £t to the CBUU-results. 





PcMS (GsV/c) 


Figure 4: Average production density of pious as a function of their CMS momentum 
for Au + Au at 1.0 GeV/A taking into account all events (0 < 6 < 16 fm); squares: 
results of the CBUU-calculation, solid line: linear £t to the CBUU-results. 
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higher than that of soft pions. 

The relation between the average prodnction density and the pion momentnm in 
the CMS of the heavy-ion reaction is depicted in fig. ^ for the reaction Au + Au 
at 1.0 GeV/A taking into acconnt all events (0 < 6 < 16 fm). Qnantitatively we 
can conclnde that pions with momenta of ~ 1.0 GeV/c are emitted on average at 
O.Spo- After prodnction of the hard pions the fireball expands to ~ O.Spo before 
most of the soft pions are emitted. Note, that both densities are below normal 
nnclear matter density snch that pions carry information essentially only on the 
low density phase of the reaction. From the linear dependence of the freeze-ont 
density on the pion momentnm we can infer that the Gonlomb potential acting on 
the pions, when leaving the reaction zone, mnst also be a linear fnnction of the 
pion momenta becanse the charge density is proportional to the baryon density. 
Indeed this behavionr can be fonnd by means of onr transport calcnlation. In the 
top part of fig. ^ the absolnte magnitnde of the average Gonlomb potential felt by 
the pion at its prodnction point is plotted versns the pion momentnm in the GMS. 
The open circles indicate the potential of the vr” and the open triangles that for 
7r~^ while the solid sqnares denote the Gonlomb potential averaged over and vr”. 
The solid line corresponds to a linear fit to the averaged valnes (sqnares) in the 
momentnm range between 0.1 and 1.0 GeV/c. We find that the average Gonlomb 
potential at freeze-ont increases linearly with the pion momentnm except for pions 
with momenta below 50 MeV/c. 

In the bottom part of fig. we investigate how the Gonlomb effects on the 
pion distribntions change as a fnnction of the centrality of the heavy-ion reaction: 
the sqnares and the circles show the average Gonlomb potential of the pions at 
freeze-ont as a fnnction of their GMS momentnm for peripheral and central colli¬ 
sions, respectively. The solid and dashed lines represent the corresponding linear 
fits to the valnes obtained from the GBUU-calcnlation. Fig. ^ clearly indicates that 
the Gonlomb effects on the final pion distribntions are also sensitive to the reaction 
geometry, i.e. pions created in peripheral collisions feel a ~ 15% stronger Gonlomb 
potential than those prodnced in central collisions. While in central collisions al¬ 
most all nncleons participate in the reaction, get compressed and collectively expand 
again, more and more nncleons do not participate in the heavy-ion reaction as the 
collision becomes more peripheral with increasing impact parameter. As a conse- 
qnence the pion freeze-ont confignration is more compact in peripheral reactions 
than in central collisions snch that the pions see a stronger Gonlomb field when 
deconpling from the baryons. 

By establishing a correlation between the momenta of pions emitted in heavy- 
ion collisions and their average prodnction times and densities we have proven that 
pions are a sensitive probe to the expansion phase of heavy-ion reactions. Obvionsly, 
by measnring the pion distribntions in heavy-ion collisions one cannot directly de- 
dnce their freeze-ont density or their freeze-ont time, bnt one can take advantage of 
the sensitivity of the pion distribntions to the Gonlomb potential dne to the presence 
of charged baryons and mesons in the reaction zone by investigating the differences 
in the observed tt’*'- and 7r“-spectra. This can be achieved by looking at the ratio of 
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Figure 5: Average Coulomb potential of pious (at their production time) as a func¬ 
tion of their CMS momentum for Au + Au at 1.0 GeV/A. Top: The circles represent 
the results obtained for tt” and the triangles those for tt’*'. The hlled squares indicate 
the Coulomb potential averaged over tt’*' and tt” while the solid line is a linear £t 
to these values for pion momenta between 0.1 and 1.0 GeV/c. Here all events are 
taken into account (0 < b < 16 fm). Bottom: The squares indicate the result for the 
Coulomb potential averaged over and tt~ for peripheral collisions {b > 6 fm) and 
the circles represent the same for central collisions (0 < 6 < 6 fm). The solid and 
the dashed line represent the corresponding linear fits to the CBUU-calculations. 
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the cross sections for 7r+- and vr -prodnction. This ratio 


R = 


) 

(j(7r+) 


( 7 ) 


is shown in £g. ^ as a fnnction of the pion kinetic energy in the CMS of the heavy- 
ion collision for the reaction Au -|- Au at 1.0 GeV/A. The top of £g. shows the 
results obtained with no Coulomb effects included while the bottom part depicts 
R obtained in a calculation including Coulomb effects. Besides the (numerical) 
statistical fluctuations the ratio R is constant (~ 1.85 ± 0.1) when neglecting the 
Coulomb forces since the pions propagate freely out of the nuclear medium. Hence 
the n~ and tt■''- spectra are not shifted relative to each other and the deviation of 
the ratio R from 1 reflects the isospin asymmetry of the Au -f Au system which 


according to reads, 


5N^ + NZ 
5^2 + NZ ’ 


( 8 ) 


where N and Z are the number of neutrons and protons, respectively. For a Au + Au 
collision eq. (|^) gives R= 1.95, which is in good agreement with our transport cal¬ 
culation without Coulomb forces. The ratio obtained when including the Coulomb 
forces (bottom part of £g. |^) shows an energy-dependent behaviour. For kinetic 
energies of ~ 50 MeV it yields a value of about 3.0 whereas with increasing pion ki¬ 
netic energy the 7r“/7r'''-ratio decreases to about 0.8 for energies of ~ 300 MeV. This 
energy-dependent behaviour of the ratio R can be attributed to the Coulomb po¬ 
tential which results in opposite forces for 7r+ and vr” mesons and which on average 
has a different strength for pions in different kinematical regions (cf. fig. ^). 


3.2 Sensitivity of the TT /tr^-ratios to the reaction geometry 

In the previous section we have shown that the variation in the Coulomb potential 
felt by pions due to their different production densities is the origin of the energy- 
dependent 7r“/7r■'■-ratio R. In this subsection we will investigate the sensitivity of 
R to the reaction geometry of the heavy-ion collision. In this respect we display 
in £g. 1^ the ratio R for the reaction Au + Au at 1.0 GeV/A as a function of the 
pion CMS kinetic energy for different CMS-angles of the outgoing pions. The solid 
line shows a fit to the total ratio i?, while the dashed, the dotted and the dash- 
dotted lines represent the resulting 'k~ ratio for Qcms = 0° ± 22.5°, 45° ± 22.5° 
and 90° ± 22.5°, respectively. We observe a large ratio for soft pions emitted in 
forward direction while the ratios for slow pions emerging at Qcms = 45° or 90° 
are signihcantly smaller. This large difference in R vanishes if the pions become 
harder. For pions with CMS kinetic energies of ~ 250 MeV the effect has almost 
disappeared. This angular dependence of the ratio R can be understood within 
the ’participant-spectator’ picture. While in central collisions almost all target and 
projectile nucleons participate in the collision and therefore form expanding hot 
nuclear matter of densities below pq, in peripheral collisions the spectator matter 
provides more compact sources for the Coulomb potential (cf. sect. 3.1). 
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Figure 6: The ratio R as a function of the CMS pion kinetic energy for Au + Au 
at 1.0 GeV/A obtained in a CBUU-calculation without Coulomb effects (top) and 
in a calculation including the Coulomb effects (bottom). The squares denote the 
corresponding computational results and the errorbars account for the statistical 
uncertainty. 
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Figure 7: The ratio R for the reaction Au + Au at 1.0 GeV/A as a function of the 
CMS pion kinetic energy for different CMS polar angles: Qcms = 0°±22.5° dashed 
line, Qcms = 45° ± 22.5° dotted line, Qcms = 90° ± 22.5° dash-dotted line. The 
solid line represents the ratio R integrated over all angles. 


The dramatic increase in the ratio R for soft pions moving in forward direction 
is due to the fact that these pions move in the same direction as the spectator matter 
and therefore stay for a long time relatively close to the spectators. This results in 
a larger Coulomb distortion of the pion-spectra than for those pions leaving the 
reaction zone in other directions. This argument is supported by £g. where in 
the top part we show the ratio R for pions moving in forward direction {Qcms = 
0° ± 22.5°) for Au + Au at 1.0 GeV/A as a function of the pion kinetic CMS- 
energy for different event classes. The dotted histogram represents the resulting 
R taking into account only central (0 < 6 < 6 fm) events, the dashed histogram 
corresponds to R obtained for peripheral events {7 < b < 16fm), while the solid 
histogram corresponds to the ratio when including all events (0 < 6 < 16 fm). 
Comparing the three histograms in the top of fig. |] we find a strong sensitivity 
of R to the impact parameter of the heavy-ion collision which is directly related 
to the impact parameter dependence of the pion Coulomb potential (cf. fig. 

For soft pions we obtain a larger value of R for peripheral events than for central 
events. The size of this effect is again reduced for harder pions. This variation of 
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Figure 8; The ratio R for the reaction Au + Au at 1.0 GeV/A as a function of 
the CMS pion kinetic energy for pions with different event classes: all events: solid 
histogram; central events (0 < 6 < 6 fm): dotted histogram; peripheral events 
{7 < b < 16fm): dashed histogram. Top: Pions moving in forward direction (0° < 
^CMS < 22.5°); middle: pions leaving the reaction zone with 45° < Qcms ^ 135°; 
bottom: all pions (0° < Qcms < 180°). 
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R with the impact parameter of the heavy-ion collision supports the argument that 
the Coulomb attraction/repulsion for pions moving in forward direction is mainly 
due to the spectator matter. The larger amount of spectator matter in peripheral 
collisions leads to larger Coulomb effects than in central collisions. This dependence 
of the Coulomb-ratio on the impact parameter of the heavy-ion reaction is found to 
be reversed for pions with kinetic energies below 200 MeV when considering pions 
leaving the reaction zone under 45° < Qcms ^ 135° as shown in the middle of £g. 

The dependence of R on the impact parameter vanishes within the statistical 
uncertainties of our transport calculation when evaluating the Coulomb ratio R for 
all pions produced in the reaction (bottom part of fig. ||). This can be interpreted 
in the same line of argument: slow pions produced in peripheral collisions moving 
perpendicular to the direction of the spectators feel weaker Coulomb forces than 
those produced in central collisions because the spectator matter moves away from 
the reaction zone and thus reduces the effective Coulomb potential ’seen’ by the 
pions. 


3.3 Comparision with experimental data 


In this subsection we compare our calculations for the inclusive pion spectra to the 


respective experimental data |W] and present 'k ~ratios for Au + Au at 1 GeV/A 
and Ni -|- Ni at 2.0 GeV/A. 

First we compare our calculation with the inclusive experimental tt■''- spectrum 
for Au -|-Au at 1.0 GeV/A for lab. angles in the range 40° < Qiah < 48° in fig. 

The full squares represent the experimental data of the KaoS collaboration as 
a function of the laboratory momentum piab whereas the histograms are the result 
of our calculation including the uncertainty due to the numerical statistics (shaded 
area). We overestimate the experimental spectrum in the range piab ~ 300-400 
MeV/c for the heavy system Au -|- Au, but reproduce the high momentum part of 
the spectrum reasonably well within the numerical statistics achieved. This general 
trend is also found in the IQMD calculations by Bass et ah [§ for this system. 
We note, that the high momentum tail of the spectrum basically stems for higher 
baryon resonances as discussed in more detail in [|n . 

In fig. 1^ we show R for Au + Au at 1 GeV/A as a function of the pion kinetic 
energy in the GMS integrated over all angles (top), for Qiab = 45° (middle) and 
©Zafe = 90° (bottom). Whereas the ratio R drops almost linearly with the pion 
energy up to about 0.25GeV when integrating over all angles, its variation with 
the pion energy is more pronounced for the two angular cuts. Here, especially the 
7T~/n^ ratio at low pion energy (< 0.1 GeV) changes signihcantly with laboratory 
angle; its drops by about 20% for Qiab ~ 45° and increases by about 20% for Qiab ~ 
90° relative to the angle averaged value. Since the relative variation in the ratio R is 
quite pronounced, it should also clearly be seen in the respective experimental pion 
ratios. For comparison we also display in £g. 0 the experimental tt / 7r’''-ratio from 
the FOPI-collaboration pl| (top: open circles) as well as the ratio obtained by the 
KaoS-collaboration for 40° < Qiab < 48° (middle: open triangels). Whereas the 


ratio resulting from the transport calculation underestimates the FOPI-data in the 
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Figure 9: The inclusive 7r+-spectrum for Au +Au at 1.0 GeV/A for lab. angles in 
the range 40° < Qiab < 48°. The full squares represent the experimental data of 
the KaoS collaboration [^] as a function of the laboratory momentum piab whereas 
the histograms are the result of our calculation including the uncertainty due to the 
numerical statistics (shaded area). 
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region of Tqms ~ 75 MeV it agrees within the statistical errors for higher kinetic 
energies. In comparison to the KaoS-data we find that the ratio is described quite 
well at Tcms = 50 MeV but is underestimated for higher pion energies. 

In fig. ^ we present the tt~/tt~^ ratio for Ni + Ni at 2.0 GeV/A for 85° < &cms < 
135°. The solid squares represent the results of our transport calculation and the 
open circles correspond to the 7r“/7r’''-ratio obtained by the FOPI-collaboration for a 
beam energy of 1.93 GeV/A [^]. Within the statistical errors of the calculation the 
data are reproduced for kinetic energies up to ~ 500 MeV. Since Goulomb effects 
decrease with the charge of the colliding system the ratio R is close to 1 when 
neglecting the Goulomb forces because N ^ Z. Due to the lower charge of the 
system the ratio shows only a very moderate enhancement for low pion energy. 


4 Summary 


In this paper we have explored the possibility to use differential 7r“/7r+ ratios as a 
clock for the expansion dynamics of a heavy-ion reaction. The analysis is performed 
within the coupled-channel (GBBU) transport approach which has proven in ref. 


nil to adequately describe differential pion spectra in the SIS energy regime. We 


have shown explicitly the correlation between the pion final production time and 
the baryon density or Goulomb potential that is seen by the pions at freeze-out. 
Furthermore, a sizeable correlation between the pion production time and the pion 
momentum is found due to the strong interaction dynamics. 

The different Goulomb potentials seen by slow/fast pions in different directions 
with respect to the beam axis and for different impact parameter of the heavy-ion 
collisions directly induce sizeable variations in the ratios especially for low 

pion momenta. These ratios can be measured directly by the FOPI and KaoS 
collaborations at GSI; detailed predictions for Au + Au at 1 GeV/A have been made 
in this context and are awaiting further experimental control. 

So far, we have performed our analysis without including any selfenergies for 
the pions while assuming the mean-field potentials to be the same for all baryons. 
From the exploratory study in ref. [23| we know, that selfenergy effects show up 


dominantly at low pion momenta in the GMS. However, since such hadronic selfen¬ 
ergies are approximately the same for tt’*' and 7r“ mesons, their effect should cancel 
out when considering the ratio of the differential spectra. Precise experimental data 
will be very helpful in clarifying this question. 
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Figure 10: The inclusive tt ~ratio R for Au +Au at 1.0 GeV/A and different 
angular cuts in the laboratory frame as a function of the pion kinetic energy in 
the CMS (squares). The errorbars indicate the numerically accuracy within the 
statistics reached. The open circles represent the experimental data of the FOPI- 
collaboration and the open triangles indicate the experimental data from the 
KaoS-collaboration pO for 40° < Qiah < 48°. 
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Figure 11: The inclusive ti~ ratio R for Ni + Ni at 2.0 GeV/A as a function 
of the pion kinetic energy in the CMS for 85° < Qcms < 135°. The solid squares 
repesent the results of the CBUU-calculation; the errorbars indicate the numerically 
accuracy within the statistics reached. The open circles represent the data obtained 
by the FOPI-collaboration for a beam energy of 1.93 GeV/A. 


20 


































[1] A. Agakichiev et al., Phys. Rev. Lett. 75 (1995) 1272. 

[2] M. A. Mazzoni et al., Nucl. Phys. A566 (1994) 95c. 

[3] G. Q. Li, C. M. Ko, G. E. Brown, Phys. Rev. Lett. 75 (1995) 4007. 

[4] QUARK MATTER 95, Nucl. Phys. A590 (1995) . 

[5] W. Gassing, W. Ehehalt and G. M. Ko, Phys. Lett. B363 (1995) 35. 

[6] Gy. Wolf, W. Gassing, U. Mosel and M. Schafer, Nucl. Phys. A517 (1990) 615. 

[7] Gy. Wolf, W. Gassing, and U. Mosel, Nucl. Phys. A552 (1993) 549. 

[8] Gh. Mntz, PhD Thesis, GSI Report 93-43 (1993). 

[9] St. A. Bass, G. Hartnack, H. Stocker and W. Greiner, Phys. Rev. C51 (1995) 
3343. 

[10] B. A. Li, Phys. Lett. B346 (1995) 5. 

[11] S. Teis, W. Gassing, M. Effenberger, A. Hombach, U. Mosel and Gy. Wolf, Z. 
Phys. A, in press; |nucl-th/9609009| . 

[12] G. F. Bertsch and S. Das Gupta, Phys. Rep. 160 (1988) 189. 

[13] W. Gassing, K. Niita and S.J. Wang, Z. Phys. A331 (1988) 439. 

[14] W. Gassing, V. Metag, U. Mosel and K. Niita, Phys. Rep 188 (1990) 363. 

[15] R. Stock, Phys. Rep. 135 (1986) 259. 

[16] K. Weber, B. Blattel, W. Gassing, H.-G. Donges, V. Koch, A. Lang, and U. 
Mosel, Nucl. Phys. A539 (1992) 713; K. Weber, B. Blattel, W. Gassing, H.-G. 
Donges, A. Lang, T. Maruyama and U. Mosel, Nucl. Phys. A552 (1993) 571 

[17] R. S. Varga, Matrix iterative analysis, Prentice Hall 1962. 

[18] J. Stoer and R. Burlirsch, Vol. 1, Springer-Verlag, Berlin, (1978). 

[19] G. Mntz et al., Z. Phys. A352, (1995), 352. 

[20] G. Mntz et ah, Proc. of the international workshop XXHI on gross properties 
of nuclei and nuclear excitations, Hirschegg, Austria, (1995), p. 345. 

[21] D. Pelte, E. Hfele, D. Best et ah, Z. Phys. A, in press. 

[22] D. Pelte et ah, Z. Phys. A, in press. 

[23] W. Ehehalt, W. Gassing, A. Engel, U. Mosel and Gy. Wolf, Phys. Lett. B298 
(1993) 31. 


21 




